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Introduction 



I pick up the discussion of quantum impurity problems where it is left in my Lecture Notes form 
the 1998 Les Houches Summer School |Q] - consult these notes for the necessary physical and technical 
background. I have so far explained how integrability provides a convenient basis of quasiparticles which 
scatter nicely among one another and at the impurity. I have showed how to use this basis to compute DC 
transport properties with a combination of thermodynamic Bethe ansatz and Landauer type arguments. 
To proceed, I now would like to discuss correlation functions. 

It is a natural idea to compute correlators using the same basis of quasiparticles, and inserting 
completude relations in the Green functions of interest. In this approach, the key ingredient is the matrix 
elements of the physical operators. Despite the "almost free" structure of the DC transport properties, 
a rather unpleasant surprise is met here however: these matrix elements are all non zero in general, that 
is, physical operators (eg the current) are able to create/destroy arbitrarily large numbers of particles 
when they act on an arbitrary state. 

These matrix elements - usually called form-factors - turn out to be a major topic |^, Q in recent 
developments (for even more recent works, see |^) about integrable quantum field theories, and are known 
for a large variety of integrable masssive quantum field theories. 

The idea to compute correlators in integrable quantum impurity problems is first to extend the form- 
factors approach to the case of theories which are massless in the bulk. The next step is to take the 
impurity (boundary) interaction into account properly, and finally, to sum over all the intermediate 
states to obtain physical correlations. The last step is the most difficult, and cannot be done analytically 
for the moment. Fortunately, and somewhat surprisingly, it turns out that the sums over intermediate 
states converge very quickly, everywhere along the trajectory between the UV and IR fixed points. As 
a result, the form-factors approach gives extremely accurate results (and with a controlled accuracy) for 
most correlators of interest. 

In some cases, the form-factors expansion is plagued with an "infrared catastrophy" due to the pro- 
liferation of soft modes in the massless scattering description. This problem can usually be controlled by 
simple means, and meaningful correlators still obtained. 

The main caveat of the form-factors approach is that it seems badly behaved when one approaches 
isotropic limits: as a result, although correlators for the anisotropic Kondo model up to g = 3/4 or so are 
now determined, the Kondo limit is still unaccessible. This is not much of a problem in applications to 
the fractional quantum Hall effect, where the interesting filling fractions are far from the isotropic limit 
of the underlying field theories. 



1 Some generalities on form-factors 

I follow here the notations and conventions of in particular section 5 of the latter notes. 
The space of states is generated by the vectors 

l"l'---'"")ai,...,a„ ==^li("l)---^l„("«)|0> (1) 

and the dual space by 

'^"'■■■''^1 («„,... = (0|Z""(a„)...^"H«i) (2) 

Here, the labels a stand for instance for solitons/antisolitons (±) and breathers in the sine-Gordon model, 
a denotes the rapidity, and the Z, are the annihilation creation operators of the Faddeev Zamolodchikov 
algebra . For a given set of rapidities and quantum numbers, the states with different orderings are 
related through S matrix elements. One usually call "in" states the states for which ai > a2 ■ ■ ■ > ctn, 
and "out" states those for which ai < q;2 . . . < a„. In or out states form a complete basis of the set of 
states, but it is convenient to keep some redundancy and consider all possible rapidity orderings, after 
division by the appropriate degeneracy factors. 
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Form-factors (more correctly, "generalized form-factors" since no order of the rapidities is prescribed) 
of an operator O in a bulk theory are defined as : 

/(ai, a„)ai,...,a,. =< Q\O[0MZlM) ■ ■ ■ <.K)|0 > (3) 

where |0 > is the ground state. We chose to take the operator at the origin, since relativistic invariance 
determines trivially the dependence of the matrix elements on the coordinates of insertion (we set h = I) 



< 



0|O(x,t)|Zt^(ai)...Ztja„)|0>=e^(^^-^*)/(ai,...,a„),,.....a„ (4) 



5* matrices have in many cases been obtained by explicit solution of properly regularized quantum 
field theories. The strategy, which was implemented quite succesfuUy in the massive Thirring model for 
instance |^], is to find out eigenstates, fill the ground state, determine the possible excitations over it, and 
finally compute their scattering, which is the result of a combination of the bare scattering and "dressing" 
effects coming from interaction with the ground state. Needless to say, the procedure is laborious, and 
it has proven much faster to obtain the S matrices in a more abstract way, using general axioms of S 
matrix theory, factorizability, and assumptions of maximal analyticity [|[ H, 

It is even harder to obtain the matrix elements from the solution of bare theories. A program to 
do this is under way ||^, 0, but very few results of interest for our problem have been obtained so far. 
Form-factors can however be determined (though not quite as easily as S matrices) using an axiomatic 
approach similar in spirit to what was done for the S matrices, and this is what we would like to discuss 
briefly here. 

One of the consequences of integrability is that multiple particle processes can be reexpressed in 
terms of two-particle ones. The basic objects in this approach are therefore those involving only a pair of 
particles. To start, consider the two particle S matrix: by relativistic invariance, it is expected to depend 
only on the Mandelstam variable 

s = {Pa{ai) + Pb{a2)Y = Ml + M| + 2M1M2 cosh(ai - cki) (5) 

The function S{s) is expected to be an analytic function of s, with a cut along the real axis running from 
s = (Ml -|- M2Y to s = 00, and another running from s ^ —00 to s — {Mi — M2)^. The exchange of in 
and out states corresponds to exchanging the upper lip of the right cut Im s = 0-1- with the lower lip of 
the right cut Im s — 0—. Crossing instead exchanges the upper lip of the right cut with the lower lip of 
the left cut. 

This plane with the cuts transforms into the "physical" strip < Im a < tt in terms of the variable 
a = ai — a2- The line Im a — Q corresponds to the right cut, while Im a — -k corresponds to the left cut. 
Upper and lower lips correspond to Re a > 0, resp. Re a < 0. 

In terms of the a variable, the standard properties of the 5* matrix are: 



Safalia) = Safalia) C 

silil{a) = Sir}{a) T (6) 



and 



• S{a) is real for a purely imaginary 
. Unitarity S^f^l (a)5,^^^g {-a) = 5l\ b% 

. Crossing 5^i^4(«) = ^^;g(*7r- a) (7) 
In addition of course, the S matrix is a solution of the Yang Baxter equation, which reads in components 
Sit. («i - "2)5,=,t ("1 - ("2 - «3) = ("2 - a3)5^^g (ai - a^)Sl\ll {a^ - a^) (8) 
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The space of states is conveniently built using the Faddeev Zamolodchikov algebra 

Z'^-{a,)Zl{a2) = S'^ll]{a,-a2)zl,{a2)Z<{a,) + 2n5Z5{c.,-a2) 



(9) 



for which several of the S matrix properties (but not those involving crossing) are just consistency 
relations. 

In terms of a, S* is a meromorphic function, whose only singularities in the physical strip are poles on 
the imaginary axis, corresponding in general to bound states. 

Consider now the two particle form-factor /(ai, 02)01,02- Up to some simple dimensionful terms 
(which are absent if the operator is scalar), this is also a function of the Mandelstam variable s. Its 
analytical properties are similar although a bit different than for 5*; in particular, the left cut does not 
exist for form-factors. 

From the second relation of (^, it follows immediately that 

/(ai,a2)Qi,a2 = S'°,^^("i - "2)/(a2,ai)4,a; (10) 
If ai > a2, the state \ai,a2 > ai,a2 stcitc, with, the Corresponding out stRte |q;2i^i ^02, ai- 

The 

exchange of the two on the other hand corresponds to going from the upper to the lower lip of the cut in 
the s-plane, ie from Im a = to Im a 27r in the a plane. It follows that 



/(a2,ai) a2,ai 



(11) 



Equations (101^) are the basic tool to determine form- factors. To see this, suppose for simplicity that 
the two particle S matrix is diagonal (actually, the following would also hold if it could be diagonalized by 
a rapidity independent change of basis: such is the case for the sine-Gordon model in the -i — sector). This 
would correspond, for instance, to the sinh-Gordon case, which we discuss in details below. Assuming 
the operator of interest is scalar, f{a) is thus found to satisfy the constraints 



/(«) = Sia)fi~a) 
fi2t7r + a) = fi-a) 



(12) 



The minimal solution of these equations is obtained by assuming that / does not have poles nor zeroes 
(except the threshold pole at a = 0) in the strip Im a G [0, 27r]. Consider now a closed contour C enclosing 
this strip. By Cauchy's theorem one has 



da 



In/™" (a) 



dz 



8i7r Jc sinlr(z - a)/2 



ln/(^) 



1 

8in 
1 

~ 8i7r 

Very often, the S matrix can be recast in the form 



dz 



In- 



smli{z-a)/2 f{z + 2iTT) 



dz 



sinh^(z - a)/2 



\nS{z) 



Sia) 



- exp 



dxg{x) sinh(a;Q!/i7r) 



from which one finds 



= est sinh - exp 



dx \ ' sin^[a;(i7r - a)/2TT] 
smhx 



(13) 



(14) 



(15) 
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In general, form-factors are expressed in terms of J""" and simple functions which encode the required 
pole structure, see below. 

Generalizing the two equations (^) we have, for n particles, 

y (oil , . . . , (Xi , tt^-f-l , . . . , O-n^ai ,...,ai,ai+i,...,a^ — f {^1 : • ■ • j : ; ■ ■ ■ ; ^n^ai ,...,bi-^-i,bi,...,an 

/ai...a„(ai + 2i7r, . . . , a„) = /a2...a„,ai(a2, ■ • ■ , ^n, ai) (16) 

The form-factors are meroniorphic functions of the rapidity differences aij in the strip < Im ay < 
27r, with two kinds of simple poles. 

One type of poles is called annihilation, or kincmatical, pole. Such poles are always expected, even if 
the theory has no bound states, when some of the rapidities in the form-factor differ by itt, corresponding 
physically to the presence of a pair particle-antiparticlc in the process. The residue is a form-factor with 
two fewer particles 

i ReSa'=a+i-^f{a' , a, ai, 012, . . . , an)aaai...a„ = [5l\ . . .51^- S'ai'"Q;^(ai ■ ■ ■ an\a)] /(ai, . . . , an)bi,...,b„, 

(17) 

where, we have defined the S matrix element 

St± («!■■■ = S'alZ («i - c.)S'4cl («2 - a) . . . (a„ - a) (18) 

This provides a resursive relation between form-factors with n + 2 and n particles, which proves crucial 
in determining the multiple particle form-factors. 

The other type of poles has a more physical origin, and corresponds to the appearance of bound 
states. In words, form-factors obey relations that mimic the bootstrap structure of the theory: if a 
particle appears as a bound state of two others, its form-factors are residues of form-factors involving 
these particles. In formulas, things are a bit more complicated. If c is a bound state of particles a and b, 
such that the S matrix has a simple pole 

5a1(«) « (19) 

then the form factor /(ai, a2, . . . , an+2)ai...a„a6 has a pole when Q!„+i — an+2 — iW^b^ with residue 

iReS Q„+i_a„ + 2=iu-j,/(ai, ^2, • ■ • , an+2)ai...a„ab = gafc/("l) Q^2, • ■ • , ^r!, ^n+l " iu'^b/'^)o.i...a„c (20) 

where the value of the n + 1*'' rapidity follows from conservation of energy and momentum at the three 



particle vertex ab — + c. Formulas ( 17|20|) are usually called LSZ reduction formulas, from the pioneering 
paper of Lehmann, Symanzik and Zimmermann 

In the simple case where there is only one type of particle (like in the sinh-Gordon model, see below), 
the general solution of these equations can always be written in the form 

/(ai, ...,«„) = K{au . . . , a„) n /™"("^ " «j) (21) 

where if is a completely symmetric, doubly periodic function, which contains all the physical poles. 
Finally, we point out that other form-factors can be obtained by crossing 

a[.....a'^ (a'l, • ■ • i^mlC'lai, ■ • ■ i ai . . ,a„ = /("l + ^l": • ■ ■ i «m + «7r, tti, . . ■ , an)a[,...Mra'^.ai....M„ (22) 



Let us now discuss indeed the sinh-Gordon model in more details - here, I follow closely the paper 
The action is 

1 



S 



dxdy 



(23) 
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Note that here the free boson is normahzed differently than in Q : this will prove more convenient below. 
Eq. (^ ) defines an integrable massive theory; the conformal weights of the perturbing operator are 
h = h = —g. The spectrum is very simple and consists of a single particle of mass M, and S matrix : 



Sia) = 



where 



tanhi (a - i^) 
tanhi (a + i^) ' 

25 



(24) 



B = 



1 



Observe the remarkable duality of the S matrix in B ^ 2 — B, i.e. in t; — s- 1/g. This duality is certainly 
not obvious at the level of the action, and is deeply non perturbative in nature. 

To determine the minimal form-factor, it is convenient to replace trigonometric functions with T 
functions using the basic identity 

(25) 



r(z)r(i-z) = ^ — 

sm TTZ 



together with the integral representation 

Inr(z) 

One finds then 



(z-l)e-* + 



dt 
T 



Sia) 



exp 



X cosha;/2 

Using a similar representation for lnsinha/2, we find finally 



(ia; cosha;(l — _B)/2 . 

sinh(ax/nr) 



/™"(a) =Afexp{8 



dx sinh(^) sinh[^^^^] sinh(2 



X sinh (.t) 

Here, we have put the normalization which is standard in the literature |T 



sin'^[a;(i7r — a)/2T:] 



(26) 
(27) 

(28) 



A/" = exp 



dx sinh(^) sinh[^^g^] sinh(f ) 
X sinh^(a;) 



Consider now the form- factors of the field $ itself. Since $, as well as the creation operators of the 
sinh-Gordon particle, are odd under the Z2 symmetry $ —^ — $, only form-factors witn an even number 
of particles are non vanishing. One finds the general formula fl^ 



/(ai, . . . , a2n+i) = M 



4sin^ 
F™™(i7r,B) 



(2n+l) 



2n+l P2n+l{Xi, 



(xi 



1<J 



where we introduced x = e" and the cr's are the basic symmetric polynomials : 

il<i2<---<ip 

with the convention cfq — 1 and Up = if p is greater than the number of variables. The P2n-i-i's are 
symmetric polynomials, which can be obtained by solving LSZ [|ll| recursion relations. The first ones 
read : 

Psixi, ... ,2:3) = 1 

P5(xi, . . . , Xs) = CT2Cr3 - C1CT5 

Prixi, ...,X7)^ O-20-30-4O-5 - (a'4a'| + cricr20-5<J6 + (^2^3 ~ c\(j2(7z) 

-22(0-1^6^7 + criCr2CT4CT7) + craCTsCTe) -f Cic|crf . 



(30) 
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with ci — 2 cos TrB/2, C2 ~ 1— cj. Observe that except for the overall normalization fi{g), these expressions 
are invariant in the duality transformation g ^. This is expected from the duality of the S matrix 
itself - as for the role of the overall normalization, it will be discussed later. 

The conventional normalization jl^ is = 1, which corresponds to chosing < 0|$(0)|a >= We 
shall make a different choice later on, when we consider the massless limit. 

The method can be generalized to models with several types of particles, like the sine-Gordon model: 
this will be discussed in the following sections. Of course, all the foregoing equations for form-factors have 
been established within the context of ordinary massive integrable field theories. The case of massless 
theories (first considered in a slightly different context in will be handled simply by taking the 

appropriate massless (ultra violet) limit in all the equations. This is somewhat safer than trying directly 
to formulate axioms for a massless theory per se, massless scattering presenting some physical ambiguities 
(see lU] for a detailed discussion of this point). 

We shall mostly consider physical properties related with the U{1) currents 9$; hence, our discussion 
will be centered on the form-factors of this operator. In the last section however, we give the example of 
a correlator involving vertex operators. 

We now discuss in more details the questions at stake in the case of the sinh-Gordon model. It has 
little physical interest in the present context, but is pedagogically quite useful 



2 Example: The sinh-Gordon model. 

2.1 Massless form-factors and the bulk current-current correlators. 

In most of the following calculations, we shall work in Euclidian space with x,y coordinates. Imaginary 
time is at first considered as running along x. The action and S matrix for the massive sinh-Gordon 
model were given in the introduction (p3|). 

Let us now try to describe the free boson theory as a massless limit of this model. First, recall the 
current correlators (we could use here the notation 0, for the chiral components, but since these get 
mixed in the presence of the boundary, we will not do so) 
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{z-z') 

< d,^z,z)d,,^z',z') - ,_ (31) 

[z — z'y 

To start, we wish to recover these correlators using form-factors. We thus take the massless limit A — > 
of the factorized scattering description of (|2^). As discussed in (l), we start by writing a = ±(A 4- (?) 
and take simultaneously A — > oo and M — > with m, m finite. In that case, the spectrum 

separates into Right and Left movers with respectively E ^ P = me^ and E = ~P — me^ . The scattering 
of R and L movers is still given by (|^ where a — > ±9. The RL and LR scattering becomes a simple 
phase, e^"^/'^. This phase will turn out to cancel out at the end of all computations, but is confusing 
to keep along. We just set it equal to unity in the following, that is we consider all L and R quantities 
as commuting. 

In this new description of the massless theory, we will need form factors in order to compute (^l|). By 
taking the massless limit of the formulas given in the introduction , it is easy to check that <& can 
alter only the right or left content of states; in other words, matrix elements of $ between states which 
have different content both in the left and right sectors vanish (that is, ^ — cj) + <f>l) 

Our conventions are conveniently summarized by giving the one particle form factor of the sinh-Gordon 
field : 

< 0\^{x,y)\9 >R ^ fi exp [me^ {x + iy)] 

<0\^{x,y)\9>L = tiexp[me^{x-iy)] , (32) 
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and we will use the obvious notation 



fiOi,. . . , 92n+l) =< O|$|01, . . . , 02„+l >R,...,B., (33) 

with the normalization of asymptotic states < 6\6' >r— 2ttS{9 — 9'). Here, / depends on the interaction 
strength g, but we do not indicate it explicitely for simplicity. We have the following properties : 

<O|$|01,...,02„+1 >J?,...,fl, - (<O|$|0i,...,02„+l >L,...,l)* 

< O|$|01,...,02„+l >fl,...,fl, - < O|$|02„+l,...,^l >L,...,L . (34) 

These form factors are expressed just as in the massive case 

nB 



r/Q a \ I 4sm— \ (2n+l) r> / N TT /mm Vl'i - I'i ; ,or\ 

\ ^mm(,*7r, r>) J I"- -I- 1"- 



where we introduced x = and the cr's are the same symmetric polynomials as before. 

We will now chose the normalization /i by demanding that the result (^l]) be recovered. Using form 
factors, this two point function expands, assuming Re z < Re z', as : 



d9l . ■ . d92n+l ___2 ^ _e, . , ^Ss^ + lN^ 



<0|a.$(z,z-)a.<i>(z',z-')|0>= - E/ (2')2.+;(2r+\)! "^'(^' 

X exp[m(z-z')(e'' + --- + e''"+'))] l/(^i,---,^2„+i)r (36) 



Now, by relativistic invariance, all the form factors depend only on differences of rapidities. Setting 
m{z — z') = e^°, (where 9o will in general be complex), one can shift all the /3's by to factor out, for 
any 2n+l particle contributions, a factor j^rpp'- Hence, the form factor expansion gives the result ( |3l|) 
provided fi is chosen such that 



oo 



^ /2„+i = 2g, (37) 



n=0 



where 



^'"+^"7 (27r)2"+i(2n + l)! + ' im, • • • , ^s^+i)! • (38) 

In practice, this sum cannot be computed analytically, but it can be easily evaluated numerically. The 
convergence is extremely fast with n, and for most practical purposes, the consideration of up to five 
particles is enough to get correct results up to 10~^. Similar convergence properties were observed in 
fisll in the massive case; note however that in this massless case, contributions with a higher number of 
particles are not damped-off by exponential mass terms. 

It must be emphasized that this result is very peculiar to the current operator. For most other chiral 
operators, the correct {z — z') dependence involves a non trivial anomalous dimension, instead of the 
naive engineering dimension. Hence, this dependence is not obtained term by term, as observed here, but 
rather once the whole series is summed up. Truncating the series to any finite n does not, in such cases, 
give reliable results all the way from short to large distances, unless some additional tricks are performed 
(see below). 

2.2 Current current correlators with a boundary 

Having fixed the form-factors normalization, let us now consider the theory with a boundary. The 
geometry of the problem is such that the boundary stands at a; = 0, and runs parallel to the y — it axis. 
The action is now : 



— / dx dy (a^$)V (9„$)V Acosh$ + A / dy cosh -$(a; = 0, y). (39) 
'^Trgj-oo J-oo L -I J-oo 2 
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This model is also integrable for any choice of A, A. The boundary dimension of the perturbing operator 
is X = —g. We can in particular take the limit where A — s- while A remains finite. It then describes a 
theory which is conformal invariant in the bulk but has a boundary interaction that breaks this invariance 
and induces a flow from Neumann boundary conditions at small A to Dirichlet boundary conditions at 
large A. As discussed in Q , the boundary interaction is characterized by an energy scale, which one can 
represent as Tb = me^^. Tg is related with the bare coupling in the action (|39| ) by A (x . In the 
following, since obviously changes of m (which is not a physical scale) can be absorbed in rapidity shifts, 
we set m — I. The effect of the boundary is then expressed by the reflection matrix : 



R{9) = tanh 



U TT 



(40) 



In the picture where imaginary time is along x, the effect of the boundary is represented by a boundary 
state. Following pq] we can represent it in terms of the boundary scattering matrix : 



\B >= cxp 



oo in 

-K{es-9)zl{e)zUo) 



|0 > . (41) 



In this formula, are again the Zamolodchikov Fateev creation operators, K is related with the reflection 
matrix by : 

i^W = i?(f -^?) =-tanh^. (42) 

We do not prove the expression (^ ) directly here. It follows however from the compatibility between 
the present computation, and the one we will do next, in a modular transformed point of view where the 
direction of imaginary time will have been switched. 

One can expand the boundary state into the convenient form : 

X 4(0i)...zl(0„)4(^i)---4(^")|O>- (43) 

Here, we have ordered the rapidities, but the result is the same as the unordered integral with the 
additional symmetry factor 1 /nl, the contributions from L and R scattering cancelling out after reordering. 
Observe now that by analyticity, the matrix elements of between the ground state and any state 
with at least one L moving particle are identically zero. More generally, the only non vanishing matrix 
elements of are those where bra and ket have the same L moving part. The same results apply by 
exchanging dz with dg and L with R moving particles. As a result one gets immediately two of the four 
current correlators : 

2g 



Z — Z' 



<dz^z,z)dz'Hz\z')> =-( 

<d,M^,z)d,,cj,{z',z')> =-(^, (44) 

which are identical with the ones without a boundary. 

The two other correlators are more difficult to get. Let us consider for instance : 

< 0|a2$(z, z)dz'^{z\ z')\B > . (45) 

The first non trivial contribution comes from the two particle term in the expansion of the boundary 
state : 

—K{eB ~e)< Q\d-z^{z, z)dz'Hz', z')zl{9)zl{e)\o >= 

X r ^K{0B - 0)e2» exp {e\z + z')\ . (46) 
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More generally, because \B > is a superposition of states with equal numbers of left and right moving 
particles, and dzcj), respectively dgcp act only on R, respectively L, particles, the expansion of (^) takes 
a very simple form : 

exp [{z + z') {e'^ +... + e^'^+O] • ■ • ' ^2n+i)p. (47) 

This correlation function depends on the product e^^ (2 + z'). It is scale invariant at the UV and IR fixed 
point. These correspond respectively to sending 6b to ^00, that is the coupling A in the action to or 
00, in other words Neumann or Dirichlet boundary conditions. In the first case, = 1, in the second, 
K = —1. Comparing with ( ^6| ) and ( |37| ) we find, as expected, that : 

< 0\d-Mz, z)d.'H^', z')\B >= ± ,_^^,., , (48) 

[z + z' )^ 

for Neumann, respectively Dirichlet boundary conditions. Although trivial, this result shows that the 
form factor expansion is well behaved, and allows us to study the correlator all the way from the UV to 
the IR fixed point when there is a boundary perturbation. In figures 1 and 2 we show the one particle 
(which is independent of B) and three particles contributions. We observe that indeed the convergence, 
by looking at the respective contributions, is very rapid. 

0.03 

0.02 
0.01 


-0.01 
-0.02 

-0.03 

2 4 6 8 10 

X 

Figure 1: One particle contribution. 

0.0002 
0.0001 


-0.0001 
-0.0002 

2 4 6 8 10 

X 

Figure 2: Three particles contribution for B — 1,0.1. 
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The only drawback of this expansion is that it is not suited for studying the correlation of two operators 
right at the boundary. Indeed in that case, Re{z + z') = 0, and the integrals in (47) do not converge. 
To solve this problem, we can introduce a modular transformed picture. We now consider the imaginary 
time as running along the y axis. Now the boundary is not represented as a state; rather, the whole 
space of states is different, since now we have only a half space to deal with. The asymptotic states are 
not pure L or R moving, but are mixtures. For instance, one particle states are : 

\\e>=\e>R+R{e)\e>L . (49) 

More generally, asymptotic states are obtained by adding to \6i,. .. ,9n >r,...,r all combinations with 
different choices of R particles transformed into L particles, via action of the boundary. Only the following 
two terms contribute : 

||01,...,0„ |01,...,0„ >R...M + ■ ■ ■ + R{Ol) ■ ■ ■ R{On)\en, ■ ■ ■ ,0i >L,...,L + .... (50) 

Although we used the same notation as previously, different things are meant by L,R. To make it clear, 
we now use the conventions : 

< 0\^{x,y)\9 >R= fiexp[me^{-y + ix)] 

< 0|$(a;, y)\9 >l= ^exp[TOe''(-y - ix)]. (51) 
To keep the notations as uniform as possible, we introduce the new coordinates : 

w{z) = iz = —y + ix, (52) 

so here R movers depend on w, L movers on w. The normalization fi is of course the same as before, and 
as before the LL and RR correlators do not depend on the boundary interaction. One finds : 

<0\d.^^w,w)d^,^w\iB')\0> =-(^^^ 

<0\d^<^>iw,w)d,^,^w',w')\0> =-j^-^, (53) 

where we used the fact that |i?(0)p 1. When compared with (p4[), these correlators have an overall 
minus sign due to the dimension h = l,h = (resp. h — 0,h = I) of the operators. 
Let us now consider : 

< 0|(9fi,$(u;,u;)a^'$(w',w')|0 > • (54) 

To compute it, we insert a complete set of states which are of the form (^9[). In the massless case however, 
since 9^,$ is a R operator, d^'^ a L operator, the only terms that contribute are in fact the ones with 
either all L or all R moving particles, as written in (po|). Thus, (|^ expands simply as : 

dei...d02n+i ,,e2„+i\2 



n=0 



^ / (27r)2"+i(2n+\)!^(^^ ~ ' ' ' ^^^'^^^ 



exp [{w ~ w'){e'' +... + e^^+O] l/(^i, • • ■ , ^2n+i)|'. (55) 



Observe the crucial minus sign when compared to ( |4^ ) . It occurs because in one geometry the correlator 
depends on z + z' , while in the other on w — w' . This now converges provided y > y' , even ii x = x' — 
ie the operators are sitting right on the boundary. Now, using the fact that from factors depend only on 
differences of rapidities, this expression can be mapped with ( |47| ) if we formally set 9 = 9' + i^, provided 
one has : 

K{0) ^R(^^^-9), (56) 

as claimed above. 
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To summarize, we can write the left right current current correlator in two possible ways. By using 
the boundary state one finds : 

/•OO 

<d-Mx,y)d,,^{x',y')>^ / dEg{E)ei,^[E{x + x')-iE{y-y')], (57) 

Jo 

(recall that x,x' < 0). One obtains g{E) simply by fixing the energy to a particular value in (^). 
When this is done, the remaining integrations occur on a finite domain for each of the individual particle 
energies since ~ there is no problem of convergence anymore. One then gets : 



OO rln E 



n=0 "^^ 



d0i... d92n 



(27r)2»+i(2n + 1)! £; - e«i - . . . - e^^- 
K{9b -e,)... K{9b - 02n)K [Ob - ln{E - e'^ - . . . - e«-)] 
|/[0i...02«,ln(i?-e«i-...-e«-)]f , (58) 



with the constraint X^il^i — E. The denominator might suggest some possible divergences; it is 
important however to realize that it vanishes if and only if the particle with rapidly 02ri+i has vanishing 
energy, in which case the form factor vanishes too. We can now shift the integrands to write equivalcntly : 



OO „o 



g{E) ^ eY^I 

^ — n ^ — 



dOi... d92n 1 



(27r)2»+i(2n + 1)! 1 - e»i - . . . - e^2„ 
X K{\n{TB/E) -Oi)... K{\n{TB/E) - 02„)i^ [HTb/E] - - e'^ - . . . - e"^")] 
X |/[0i...02«,ln(l-e'^ -...-e^-)]f , (59) 

where the constraint < 1 is implied, we used the fact that form-factors depend only on rapidity 

differences, and Tb = e.^" ■ 

By using the dual picture, one finds, instead of (|57|) : 

(60) 



< 9j$(x, y)d,.^{x\ y') / dET{E) exp [-iE{x + x') - E{y - y')] , 

Jo 



where : 

OO „o 



^/ 7-^\ \ ^ f dOi . . . d62n 

^ ^ = " ^J_^ (2^)2"+i(2n+l)!r 



1 

- . . . - e« 



X R{ei - HTb/E)) . . . R{02n - HTb/E))R [In (l - e^^ - ... - e''^") - ln(rB/£;)] 

X |/[0i...^^2„,ln(l-e''i -...-e«-)]f , (61) 

where in (^9|) and (^) the constraint Y^i=i < 1 is implied. The two expressions are in correspondence 
by the simple analytic continuation : 

g[E)^iT{iE). (62) 



3 The sine-Gordon Model. 

In this section we follow the same line of thought for the sinc-Gordon model. This is the massive 
deformation of the free boson which preserves integrability with either boundary interactions used in the 
fractional quantum Hall problem and the anisotropic Kondo model [|]. Thus the form factors of the 
sine-Gordon model in the massless limit will be the quantities we need. The solitons/anti-solitons and 
breathers quasi-excitations make the problem more complicated but the results presented before hold 
with the addition of a few indices (and rather more complicated form factors). 
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We start as before with the massive sine-Gordon model, whose action reads 



1 



S=-^l dxdy {d^^f + {dy'^f + K cos . (63) 



Notice again that we have used a different normahzation than in this wiU avoid carrying factors of 
2tt all along the following paragraphs. Recall that g = P the usual sine-Gordon parameter. 

The form factors approach is formally the same, albeit more complicated because the particle content 
is much richer, and depends on g. For 1/2 < 5 < 1, only solitons/anti-solitons appear in the spectrum 
of the theory. This is the so called repulsive case, with g = 1/2 the free fermion point (giving rise, in 
quantum impurity problems, to so called "Toulouse limits". When < 5 < 1/2, the particle content 
is enriched by \\/g — 2] bound states, called breathers. In the following we will denote by the indices 
a = ± the solitons and anti-solitons, and a = 1,2, [1/g — 2] the breathers. The solitons form factors 
in the massive case were written by Smirnov ||] and we obtain the massless form factors by taking the 
appropriate limit of the massive ones. Only right and left moving form factors survive in this limit, as in 
the sinh-Gordon case. Moreover, the symmetry of the action dictates that only form factors with total 
topological charge zero are non-zero for the current operator. As an example, the soliton/anti-soliton 
form factor is given by : 

<O|^a,$(z,z)|0i,02 = a' ^imd e'-'^+'-y^ 

ZTT 

^^^'-^^^ exp[m(e^^+e«^)z], (64) 



cosh ^^2^(6*1- 6*2 +i7r 



with a + a' — and a — ± stands for soliton (resp. antisoliton). From |^ one has 

g.^2 x{^^-e) gj^j^ (l-2g)x 

= csinh 1 exp I / . ^ . ^ ''''f^ rfx I , (65) 

2 \ ./n 2: smh 2(l-g) x cosh f ' 

(this is essentially the minimum form-factor discussed in the introduction) with the constant c given by : 

4(1 /l r ■""'■§ ™''T(T?Sf 



and d by : 

c^=^^, (67) 

ZTTC g 

The normalization constant /i can be determined from first principles. Indeed, the operator being 
related with the U{1) charge, we need that 



<0i|/ —a,$|02 >f =2^(5(01-02), (68) 

J -00 



using the fact that a soliton for the bulk theory (|63| ) obeys ^(cx)) — $(—00) = 27r. On the other hand, 
using the dependence of the form-factor on spatial coordinates, the left hand side is 

I < O,\^dx<!>\02 >^ / da;e""(^'^-^'')- = ^r< dil^d^m >+ m - ^2), (69) 

ZTT J me" ZTT 

Comparing the two and using crossing leads to the identity 

indCiin) = 1 (70) 
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or 



/i = 2tt- 



(71) 



1 — g cd 

The other soHton-antisoUton form factors follow from the sophisticated analysis of jsj . Their expression 
simplifies in the case g = i, t an integer. This is the physically relevant case for the v — j fractional 
quantum Hall effect. One then finds : 

1 



<o|— 5,$(z,z)|0i,...,02„ >^:::^ 



,+...+ 



sinh 



9p - in) 



Kj 

71 2n 

Jl Jl sinh~^(t - 1)(6'5 - 6lp) detiJ. 



=1 g=n+l 

The matrix H is obtained as follows. First introduce the function 

t-2 



^(a) =2*-2[]sinh- 
J=l 



a — I 



TTj vr 

h i — 

t-1 4 



One then defines the matrix elements as 

^0 k—2n 



H.. 



1 

2i7r 



-2i7r 



da Jl V(a - Pk) exp [(n - 2j - l)a + (n - 2i)(t - l)a] 



(72) 



(73) 



(74) 



fe=i 



where i, j run over 1, ... ,n — 1. It is not difficult to convince oneself that this produce a symmetric 
polynomial of the right degree. Although cumbersome, it is an easy task to extract these determinants, 
as examples we find for g — 1/3 : 



a^ie'"), n = 2, 



detH = exp{~^e^) cri(e'''')o-3(e^''), n = 3, 



(75) 



up to irrelevant phases and with the Ug's defined previously. Having these expression we can get all form 
factors using the axiomatics sketched in the introduction. For example, the solitons form factors with 
different positions of the indices we use the symmetry property (|l6|): 



■5 P^, (7^+1, t7^ 



)ai,...,a^,ai^i,.. . ,a„ 



)ai,...,6i + i,bi,...,a„"5'ai,ai+i (^i ^i+l)- 



(76) 



Here again, we omit the distinction between left and right moving form factor, they are simply related 
by complex conjugation. At the points g — l/t the soliton S matrix used in the last expression is 
reflectionless and basically just permutes the rapidities up to a phase. When there are breathers, the 



soliton S matrix has poles corresponding to the bound states at the points 9 = iu'^i^ = in 



■nrg 



m for 



the m'th breather. In view of the last relation, this induces poles in the form factors. We obtain the 
breather form factors from these poles : 



zrese„_j- 



= r„,(-l)^™(5„„_,+,„ 

X /(gi, + — -l ^^^^_^^ )a,,...,a^.,,m, (77) 



and is given by the residue at iu^^ 



b (1-9) 



C + + 



mmg 



l^gj {l-g) 



9 ■ 1 - .9 
■ smTT 



1/2 



(78) 
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Having these relations, we posess all ingredients to compute all form factors for g — l/t. Using them 
for the computation of the current correlations is then merely an extension of the previous results for 
sinh-Gordon with indices. The normalisation of the form factors, fi, should also ensure that ( ^4|) is 
reproduced. This is fixed by introducing a complete basis of states : 

and computing the correlations exactly like in the sinh-Gordon case. 

Keeping a finite number of form-factors and demanding that the two point function is properly 
normalized gives rise to values of /i which are slightly different from (|7l|). How different is a good 
measure of the convergence of the expansion, and the validity of the truncation. For g — 1/3, the one 
breather and 2 solitons form factors normalise to ^ — 3.14 which is very close to the exact n. Similarly for 
g = 1/4 we found from the contributions up to two solitons that /i = 2.05 to compare with 2.094 = 27r/3. 

Calculations in the presence of an integrable boundary interaction are also done like in the sinh-Gordon 
case. The boundary state is now given by : 

■/o<ei<...<9„ Ztt 

X Zl-(0i)---4""(^n)4'n^?l)-4''"(^n), (80) 

with an implicit sum on the indices. The matrix i^T"'' is related to the boundary R matrix in the following 
way : 

K^'iO) = (z| - e) . (81) 

The b means that we take the conjugate of the indices ie. ± ^ =p and m m. 

From the previous expressions, we can compute de current-current correlation function in the presence 
of a boundary for 5 = l/t. The results we will get depends on the boundary interaction, in the next 
subsection we present some results for the boundary sine-Gordon model, which is of relevance to tunneling 
experiments in fractional quantum Hall devices jl[ ^ ^ |l9j . 



4 Conductance in the fractional quantum Hall effect 
4.1 General remarks. 



The boundary sine-Gordon action is 



S 



1 



dx 



dy 



(a.$)2 + (a,cf>)2 



+ A 



dy cos -<^{x = 0,y). 



(82) 



The reflection matrices have been worked out in [|6|, For generic values of the coupling g, the 
amplitude for the processes H > + and > — is R±{0 — 9b), and for the processes H > — and > -|- 



it is R^{e - 
not raised): 



0b) with (the notation is slightly changed with respect to M, where outgoing indices were 



R^ie) 
Riie) 



(i-g)f , , 

= e-^R{e) 

(g-l)e 

= ie~^R{e) 



(83) 



where the function R reads 

m 



2 cos 



g)e _ 



n 



1=0 Y,(-e) 



r(j+i^-^^)r(i+('+i)il 



-a) i(l-3)9 ^ 



r ( i +0+1/2) ^ - ii^) r ( I +(i+i/2) - i<i^ 



(84) 
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In (|83|), our conventions are such that in the UV hniit {9b — oo) the scattering is totaUy ofF-diagonal so 
a soUton bounces back as an anti-sohton, in agreement with classical limit results for Neumann boundary 
conditions. A useful integral representation of R is given by : 



o«7 



2 cosh 



dy . 2(1-5)% smhi^)y 

— sm — 

2?/ ffTT sinh 2y cosh ■ 



R(9) = f J- exp i / — sin ^ — - . (85) 

- ra_-g)e_.^] ^\ 2y gn sinh 27/ cosh ' 

2g 



Recall that the spectrum is made of one breather and the pair soliton antisoliton in the whole domain 
1/3 < g < 1/2. More breathers appear for g < 1/3 (and the reflection matrix of the 1- breather is always 
the same as in the sinh-Gordon case.) There are no breathers for 5 > 1/2. 

The physical quantity of interest in this case corresponds to the AC conductance at vanishing tem- 
perature in the edge states tunneling problem: this problem has been described in details in |l[] , to which 
we refer the reader. 

A standard way of representing the conductance is through the Kubo formula ||2l| : 

-1 /'L /"OO 

G{lom) = -^:^^j^ j ^dxdx' j dye'^^^y<j{x,y)j{x',0)>, (86) 

where com is a Matsubara frequency, y is imaginary time, y ~ it. One gets back to real physical frequencies 
by letting ojm = —iuj. In (|86| ) , j is the physical current in the unfolded system |l|. Without impurity, 
the AC conductance of the Luttinger liquid is frequency independent, G — g. When adding the impurity, 
it becomes C? = ^ + AG. After some simple manipulations using the folding §, one finds : 

/.oo 



AG(a;A/) = z / / dye 

SttlvmL'^ J-l J-oc 



X [<dM^,y)dz'Hx',0)> + <d-M^,y)d,,<S>ix',0)>], (87) 

where z = x + iy. The strategy is simply to evaluate the current-current correlator using form-factors, 
and extract the conductance from (p7|). 

We will do so for special values of g, but first, we can extract some general features of the UV 
and IR expansions easily. To do so, consider the soliton antisolitons reflection matrix. Evaluating the 
integral in ^ by the residues method leads to a double expansion of the elements in powers of 
exp(0) and exp(i — 1)6* . This leads for the conductance to a double power series in {uj/Tb)^^^^^^ and 
{uj/TbY in the IR, {Tb/ojY~'^^ and [TB/i^Y in the UV. Breathers do not change this result. For instance 
for the 1-breather, since the reflection matrix is the same as in the sinh-Gordon case, and therefore g 
independent, the contributions expand as a series in [uo/TbY ™ the IR, {Tb/ojY in the UV. This holds 
for any coupling g. Therefore, as first argued by Guinea et al. at low frequency, the conductance 

goes as oj 2 for 5 < 1/2, for 5 > 1/2. The power would seem to indicate that there should 

be a term in the DC conductance, but this is not correct because only the modulus square of i?^ 
contribute to the DC conductance, and this expands only as powers of exp(i — l)/3 . 

The presence of analytic terms in the IR is a straightforward consequence of the fact that IR per- 
turbation theory involves an infinity of counter-terms, in particular polynomials in derivatives of $ p2| . 
More surprising maybe is the fact that we find analytical terms in the UV. This requires some discussion. 
The UV terms follow from the short distance behaviour of the correlation function of the current. For 
any operator O we could write formally, 

00 „ ^ ^ 

< 0{x', y')0{x, y) >x= ^(A)^" d\...dn< 0{x' , y')0{x, y) cos -$(1) ... cos -$(n) >a=o ■ (88) 

From (p8|), one would naively expect that the two point function of the current expands as a power series 
in A^, which would lead to a power series in (tj/Tg)^^"^. This is incorrect however because, even if 
integrals are convergent at short distance for g < 1/2, they are always divergent at large distances. It is 
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known that these IR divergences give precisely rise to non analyticity in the couphng constant A. One 
usually writes : 

< 0{x',y')0{x, y) >= ^ C^oi^' - x,y' ~ y)0,{x, y), (89) 

i 

where Oi are a complete set of local operators in the theory and the C's are structure functions. These, 
being local quantities, have analytic behaviour in A. However, < Oi{x,y) > being non local is in general 
non analytic - actually, on dimensional grounds, 

< 0,(a;,y) >cx A^/(i-f) cxTI', (90) 

where A = h + h is the (bulk) dimension of the field Oi . If we computed the conductance perturbatively 
using Matsubara formula, we would use (|8|) with O the electrical current operator. The case Oi the 
identity operator gives rise to an analytical expression in A, but eg the case Oi = dzdg^ gives \'^/(^^~3) 
times an analytical expression in A (its mean value can be non zero because there is a boundary). More 
generally, since the only operators Oi appearing in the case of the current are polynomials in derivatives 
of <i>, all with integer dimensions, we expect that the two point function of the current will expand as a 
double series of the form A^"A^™/^^~^'', ie going back to Tb variable , that the conductance will expand 
as a double series of the form (rB/w)2"(i-9) (rB/w)2™, in agreement with the form factors result. 



4.2 The free case. 

In the case g — 1/2 one has simply 



Rt{0)^Q{e) =-^. (91) 



In that case, only the soliton-antisoliton form factor is non zero, f {91,02) = i/ie^i/2ge2/2^ ^-^^^ ^^^^ 
normalization fi = 2tt and we have set m = 1. Hence, J^{lu) is readily evaluated 



oo roo 



T{oj) = / / deid926{e''' + e'^ ~ io)[Q{0i)Q{92) - P((?l)P(^?2)]e"^e*^ (92) 



oo — oo 



so 

T{uj) ^uj f' dx- ^^^J'^^^^ ^ (93) 







from which it follows that 



Thus we find 



2 V 

This is in agreement with the solution of [l^ 



AG{u;)^--^t^n-\iu/TB). (94) 
G{uj) = i f 1 - ^ tan-i {uj/Tb)] ■ (95) 



4.3 G{uj) at g = 1/3. 

The conductance for g — l/S has a direct application to the quantum Hall effect. Comparing with the 
free case, previously treated, we now have a breather in the spectrum and non zero form factors for all 
number of rapidities. Still the convergence is such that evaluating the first few form factors give results 
to a very good accuracy, independently of the regime, UV or IR, in which we make the computation. 

In this case, the first few non zero form factors are /i, /±.ip, /±.=p.i, /i,i,i, etc... Here the subscript 
"1" denotes the breather. The first step is to compute the normalisation in order to satisfy (^l|). When 
computing this normalisation, we find that the first two form factors account for the whole result to more 
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than one percent accuracy. Then including the 1 breather-2 solitons form factor is sufficient to get the 
result to a very good accuracy is negligeable) . Actually one observes that the speed of convergence 

of the form factor expansion varies geometrically with the number of solitons (counting the breathers as 
2 solitons). 

In order to get the conductance we need the reflection matrices, they were given in previous expressions 
and reduce to a simpler form for this value of g : 



m 



1 



r(3/8-f )r(5/8 + |^) 



2cosh(e'- f ) r(5/ 



:)r(3/8 + f) 



and the breather reflection matrix is 



Rl{e) = tanh 



in 



(96) 



(97) 



From the pole of the 2 solitons form factor, the one breather form factor is found using ( [TT] ) and its 
contribution to the conductance is : 



AGicjy-' = ~fi^ 7^e tanh 

8 



log( 



i7r/4 



(98) 



here fi = n is fixed by (|71|) and d 
computed similarly, we find : 



0.1414.... The contribution from the two solitons form factors is 



(2) 



cosh^(e-log(l-e»)) 



(99) 



where (^{9) is the function defined in (|6^). We can similarly write the following contribution, and we 
find that these last two expressions are sufficient for any reasonable purpose, they give the frequency 
dependent conductance to more than one percent accuracy. We give the full function G{lu) in figure 3. 




0.001 0.01 



10 100 



Figure 3: Frequency dependent conductance at T=0. 



Observe that in the UV and in the IR we obtain the lo dependance discussed previously, even with the 
truncation to a few form-factors. The form-factors expansion is indeed very different from the perturbative 
expansion in powers of the coupling constant in the UV, or in powers of the inverse coupling constant 
in the IR. Each form-factor contribution has by itself the same analytical structure as the whole sum; 
contributions with higher number of particles simply determine coefficients to a greater accuracy. 
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5 Anisotropic Kondo model and dissipative quantum mechan- 
ics. 



It turns out that the anisotropic Kondo model is related with the famous two state problem of dissipative 
quantum mechanics [E3| . The bosonised form of the hamiltonian is : 



2 



dx 



Sirg 



■i$(0)/2 



(100) 



where S are Pauli matrices. As for the quantum Hall problem, we keep using as a basis the massless 
excitations of the sine-Gordon model; however the boundary interation is different: this will result in 
different reflection matrices. 



5.1 A computation of C(t). 

We first work in imaginary time. We consider therefore the anisotropic Kondo problem at temperature 
T. Let us consider the quantity X{y) =< [S^{y) ~ 5*^(0)]^ >. On the one hand, using that S*^ = ±1, it 
reads 2[1 — C'(y)], where C(y) is the usual spin correlation 

Chi) ^\{< S\v)S\^) > + < S\^)S\y) >] . (101) 

On the other hand, we can write a perturbative expansion for X(y) by expanding evolution operators in 
powers of the coupling constant A. At every order, we get ordered monomials which are are a product of 
a monomial in 5^ and vertex operators of charge ±1/2. We must then evaluate S'^(y) — 5^(0) for each 
such term, trace over the two possible spin states, and average over the quantum field. Since we deal 
with spin 1/2, terms S'^ and S~ must alternate, and there must be an overall equal number of S""*" and 
S'~, and an equal number of 1/2 and —1/2 electric charges. 

Now, since each S~^{y) comes with a e~'*'^'/^ and each S~{y) comes with a e'*'^'/^, S^{y) = 5^(0) if 
there is a vanishing electric charge inserted between and y, and S^{y) = —5^(0) if the charge inserted 
between and y is non zero (and then it has to be ±1/2). Therefore, we can write the perturbation 
expansion of X{y) in such a way that the spin contributions all disappear: 

I °° 2n y y y l/T „l/T 

^{y) = vX!-^^" Y dyi dy2... dyp dyp+i . . . dy2n 

n=0 alternating ei=± p=0-'° "'wi "'fp-i Jy2,^-l 

X 4(ei + ... + ep)2{Je-"i'^(2'i)/^..e-"^'"^(''^")/2^^, (102) 

where Z is the partition function, the factor 4 occurs because of the normalization 5^ — ±1, for every 
configuration of e's, only one value of 5^(0) gives a non vanishing contribution. Here, the label N indicates 
correlation functions for the free boson evaluated with Neumann boundary conditions (the conditions as 
A^O). 

On the other hand, let us consider the correlator 

< 5,$(x, y)$(0, y') >n^ "85^—7^ ^, (103) 

± (y - y'y 

which goes to —8g7r6{y — y') as a; 0. We have then, by Wick's theorem, 

< e-"i*(''i)/2 . . . e-^^^"*(''^")/25,$(x, y) >n = 

8*5 I >Jq ^, , ^. ) (e-i*(''i)/2...e'^-*(''-/2))^, (104) 
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and therefore 

C2n \ / 2n \ 
H ^^^^ ^ E ^^^^-7 ^ ( . . . , (105) 



where contractfons between the dots : are discarded. In (105), contractions between the dots would lead 
to a term factored out as the product of the two point function of dx^ and the 2?! point function of 
vertex operators, both evaluated with N boundary conditions. Now, we are going to be interested in the 
X — > limit where, with N boundary conditions, dx'^ vanishes. As a result we can actually forget the 
subtraction in 105, and write simply obtain 

X{y, A) = -j-^ lim f dy'dy" < dx^x, y')dxHx, y") >a, (106) 

where the label A designates the correlator evaluated at coupling A, N corresponding to A = . Hence, 
we can get C{y) from the current current correlator. The latter can then be obtained using form factors 
along the above lines. The only difference is the boundary matrix. If we restrict to the repulsive regime 
where the bulk spectrum contains only a soliton and an antisoliton, one has : 

i?^ = tanh (^^ - ^) , i?±=0. (107) 

Here again our conventions are such that a soliton bounces back as an antisoliton, in agreement with the 
UV and the IR limit that have Neumann boundary conditions. In the attractive regime we need to add 
the breathers with : 

tanh(; 



_ ^2 4(1/9-1) 

™ tanhl 

Writing, as in (|5^) 



^ ^ -AiniZ±r_ (108) 



< 5.-$(x, y')dMx, y") >A= / Q{E, (3b) exp [2Ex - iE{y' - y")] , 



(109) 



we have that 



lim^_o < dx^{x,y')dx'^{x,y") >x = 
dE [g{E, (3b) - g{E, -oo)] exp[-tE{y' - y")] + c.c, (110) 

where the < dz^dz^ > part and its complex conjugate (which are A independent) have been evaluated 
by requiring that the correlator vanishes as A ^ due to A'' boundary conditions. Hence, using the fact 
that Q is real, 

X{y) = ^ [GiE, Ob) g{E, -oo)] ,in\Ey/2). (Ill) 

Therefore, if we write : 

POO 

C{y)-1= A{ujM)cos{ujMy) duM, (112) 



where wa/ is a Matsubara frequency, we have : 

1 1 

(257r)2 w|f 

An observation is now in order. From the foregoing results we see that 



^^^m) = TTT^— [G{^m,0b) - GiiOM,-^)] ■ (113) 



<S^{0)S^{y)>-l = ^—^ —[g{E,eB)^g{E,^^)]cosEy. 



(114) 
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On the other hand, consider the expression 



< 



dx' 



dx" [< d,<fix',y)d,^x",0) >x-< a,$(x',2/)a,$(x",0) 



(115) 



By using the same representation (109), this is 

dx' / dx" / dE [g{E, 0b) - g{E, ~oo)] exp [E{x' + x") - iEy] + cc, 
-oo <y —Qo Jo 

which coincides with ( |114| ) after performing the integrations. We conclude that 

< s%o)s%y) > -I =< JMMv) > - < JMMv) >N, 

where we defined 

1 f° 

Jx^T, — / dx^{x,y)dx. 

2.971- J-oo 

We find also by the same manipulations that 

< S\^)S\y) > -I =< Jy{0)Jy{y) >-< JyiO)Jyiy) >n, 

where 



1 



dy(j){x,y)dx. 



We now continue to real frequencies to find the response function : 

x"M^i / dte^-'{[S%t),S%0)]), 



to find : 



2 

1 1 



■Im [g{-iLU, 9b) ~ g{-iuj, -oo)] 



(116) 
(117) 

(118) 
(119) 

(120) 
(121) 



As a first example, let us consider the so called Toulouse limit or free fcrmion case. Then the only 
contribution comes from the soliton antisoliton form factors, which as discussed above is f {01,62) ~ 
iiie'^^l'^eP-^l'^ . Hence, 



X"(u;) = ^Re 



00 /'OO 



^611^6*2 



1 



00 ./ — 00 



that is 



and 



X"(c^) = 4^1mf Tdx 



(e^i + ie^s)(e^2 + ie^^) e^i + e®^ 
1 



5(e^i + e^^ - cj). 



4 Tn 



-Im- 



1 



(a; + iTB){y> — a: + iTb) ) ' 
w + iTb ' 



(122) 



(123) 



TX^ UJ CiJ + HTb 



In 



1 4r| 



1 



■In 



1 



■ tan 



Tb 



(124) 



In general, observe that, since the reflection matrix for solitons and antisolitons expands as a series 



„/3 x"(") 



, will, for any coupling, expand as a series of the form (w/Tb)^" in the IR. In particular, this 



leads to a behaviour C[t) (x >> 1 for any g. In the UV, one has to split integrals in two pieces. Since 
the soliton-soliton form factors expansion involves powers of exp(i — 1)9, ^ ^"-^ expands as a double series 
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in {Tb/ujY^'^^ and {^TbI^Y ™ the UV. Hence at short times, C(t) — 1 cx This is in agreement 

with the quahtative analysis of 

Results for 5 7^ 1/2 are more involved because there are non zero form factors at all levels. Still, when 
working out the first few form factors we observe a very rapid convergence with the number of rapidities 
and again we can give precise results for different values of g. As an example, let us show the results for 
5 = 1/3. 

The computation for g = 1/3 is very similar to the previous conductance computations. The boundary 
matrices are much more simpler though. In this case we have : 



i?^=tanh(l-^), i?|=0. 



(125) 



and 



tanh(T 



(126) 



tanh(| + f ) 

Then, as was found for the conductance, we find that the first two contributions are sufficient for most 
purposes, they are given by : 



Ax"(a;) 



2j2 



■Re 



tanh 



-) 



tanh 



(127) 



and 



Ax"(^) 



fZ^idV 1 



Re 







IC(^^-log(l-e^))l%. 



,_oo cosh"((?-log(l-e«)) 
[R+_{e + \og[uj/TB))Rt{\og[[l - e%uj/TB]) - l] 



(128) 



Again these two expressions are sufficient to get a very precise result. Similar computations give rise to 
the results in figure 4 where we plotted S'(w) = \"{u))/uj for the values, g = 3/5, 1/2, 1/3, 1/4. 




Figure 4: Spectral function for Tb = 0.1. 



When making these calculations we have to be careful about which terms are needed for a good 
convergence. Our observation is that keeping the form factors up to two rapidities give very good results 
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precise to 1%. It is possible to go further and get a better precision if needed. The last statements are 
true for g e [0.6, 0.2] and we believe even further (we can get rough bounds on the higher contributions 
and have an idea of the precision). Still, for the moment the isotropic Kondo point is difficult to treat. 

To a good approximation, S{ijj) can be represented for (7 < i by a Lorentzian shape built around the 
four poles of the one-breather term: 



est 



(129) 



where the poles are located at w 
that 



±51 + ±ir. Here, it follows from the one breather reflection matrix 



f 



= cot ■ 



Tig 



(130) 



and = Tb sin 



■Kg 



2(1-5)' 

The approximation improves as g gets smaller. 
According to (129), the pic in S{lo) disappears at g = i, where £7 = F. The poles meanwhile reach 
the imaginary axis at g = ^ where p = 0. These features are shared by the exact solution, as far as can 
be established from numerical study of the form- factors series: the pic in S{ijj) disappears at g = ^, while 
the tail oscillations of C{t) disappear at 5 = i. Depending on one's definition of the "transition" from 
coherent to incoherent regime, the latter therefore occurs at g = ^ or g = i. This feature was somewhat 
unexpected from previous approximate studies, but is fully confirmed by recent numerical work j2^, [26) . 
In any case, there is no real transtion per se, and it is possible that the study of different quantities (for 
instance multi point correlations) would exhibit a qualitative change of behaviour at still another value 
of g. 



5.2 Shiba's Relation. 



Up untill now, we showed results for certain values of g more or less limited by our ability (or tenacity) to 
write the form factors corresponding to that value of the anisotropy, and make them converge. It is not 
impossible to find general relations though; for example the behaviours in the UV and the IR in different 
models were infered in all generality. 

Here we present a generalisation of Shiba's relation p7 [which was proven for the Anderson model 
and generalised to Luttinger liquids by Sassetti and Weiss 1 28| . The relation states that : 



lim 



X"(t^) 



(131) 



with xo the static succeptibility. If we look at the quantity 

dOi . . . d02n-l 



g{E) = EJ2 [ 

*i — n — 



n=0 ' 
aifci 



(27r)2»(2n)! 1 - e^'i - 

if"i''i(ln(rB/S) -Oi)... K"^-^^^ 
j^a2^b2^ [ln(rB/S) - In (1 - e^i - 



HHTb/e) 



(132) 



insert it in the expression for x"i^) ^-iid expand it around E' ~ we find that the contributions from 
the K matrices all cancel and only a constant is left (we have to take into account the fact that the 
soliton/anti-solitons K matrices always appear in pair). Then comparing this with the UV normalisation 
we find that : 

' (133) 



lim 

0)^0 



22 



The total succeptibility is x = x' + ix" and the static succeptibihty xo which is the zero frequency hmit 
of x' can also be infered from t he p revious expressions for the spin-spin correlation. We just need to take 
the real part when continuing (113) to real frequencies, which leads to : 



Xo 



TT 



B- 



(134) 



Finally, in order to make contact with the usual form of Shiba's relation, we need to renormalise the spins 



to 1/2 and use the usual normalization for Fourier transforms, which leads indeed to (131). 



6 Friedel oscillations: correlations involving Vertex operators 

As exemplified in the previous sections, the method works naively indeed for currents, ie for operators with 
no anomalous dimension (calculations for the stress energy tensor for instance, would be very similar). 
Many physical properties are however described by more complicated operators, that is operators which 
have a non trivial anomalous dimension. As an example, I would like to discuss here the equivalent of 
Friedel oscillations in Luttinger liquids: more precisely, the 2kp part of the charge density profile in a 
one dimensional Luttinger liquid away from an impurity. This is a problem which has attracted a fair 
amount of interest recently ||2^ , |3^ . 

We start with the bosonised form of the model. The hamiltonian is: 



1 1 

H^- dx [87:gU^ + —{d,(bf] 

2 7_oo J57rg 



A cos 0(0), (135) 



where we have set vp = 9- Then for the Friedel oscillations, the charge density operator is just : 

p{x) — Pq + 28x4' H cos[2kpx + 0(x)]. (136) 

TT 

with po = the background charge. 

We decompose this system into even and odd basis (this is was explained already in |^) by writing 
4> — 4>L + (f'R and setting : 

^''{x + t) = ^ [(I)l{x, t) + (Pr{-x, t)] 

ip°{x + ^) = ^ [4>L{x,t) - M-x,t)] (137) 

Observe that these two field are left movers. We now fold the system by setting (prefactors here are 
slightly different from ([|l|) because of the normalizations in the hamiltonian) : 

y/2Lp''{x + t), a; < 0^ = V2(?!)'=(-x + t), a; < 
41 = y/2if''{x + t), X <Q (j)% = -^/2(j)''{-x + t),x <0 (138) 

and introduce new fields 4'^'° = 4'l° + 4'r° ■ The density oscillations now read : 

m_J!A = eos(2fc^x + ,^)(cos ^) (cos ^), (139) 
Po 2 2 

with rjp the additional phase shift coming from the unitary transformation to eliminate the forward 
scattering term |^. 4° is the odd field with Dirichlct boundary conditions, at the origin (/>°(0) = 
leading to [|lj : 

■ 1 \ 



(cos^l-) cx ( - ) , (140) 



23 



and the (f)'^ part is computed with the hamiltonian : 



dx [STTgW^ + - — {d^cj}")'^] + Acos ^ 
8ng 2 



(0) 



On general grounds, we expect the scaUng form 



9/2 



(141) 



(142) 



where -F is a scahng function to be determined. Note that even the smah x behaviour of this function 
was not known in general. 

To compute the correlation functions , we use the formalism introduced in the previous section (with 
(j)^ = $). The massless scattering description and the boundary state are the same; only the operator 
which is studied changes. 

The one point function of interest is (0| cos ■§ |-B). We thus need the matrix elements of the operator 
cos ^ in the quasiparticle basis: these still follow easily from the massive sine-Gordon form- factors ||], 
like the form- factors of 5$ in the previous sections. Difficulties however arise when one considers the 
contribution of each form-factor to the one-point function: the rapidity integrals turn out to be all IR 
divergent! This was not the case for the current operator, whose form factor has the naive engineering 
dimension of an energy, leading to convergent integrals: some sort of additional regularization is thus 
needed here. 

To explain the strategy, consider first the case g = 1/2. Here, the friedel oscillations are simply 
related to the spin one point function in an Ising model with boundary magnetic field. By using the same 
approach as the one outlined before, one finds the following form-factors expansion : 



00 „oo " 

,1=0 



— — tanh — 

2tt 2 



Jl f tanh- 



(143) 



The integrals are all divergent at low energies, when 9i —00 and the integrand goes to a constant. Let 
us then introduce an IR cut-off ( we chose 9 > 9,ni7i and set A = e^""") and take the log of the previous 
expressions (a similar method has been used in ^ to study the UV limit of massive correlators. See also 
||37||,|3^). Ordering this log by increasing number of integrations, one can show that each term diverges 
as In A. Moreover, since the divergence occurs at very low energy, where the tanh goes to unity, the 
amplitudes of these InA do not depend on 9b (for 9b ^ —00), ie on the boundary coupling. It is then 
easy to get rid of the cut-off: we simply substract the log of the IR spin function, ie we substract the 
same formal expression with 9b — 00. The first two terms of the resulting expression read : 



In 



{a{x)) 



IR 



du 
27ru 



Tb - u 
Tb + u 



1 



-j-j- dUj ^_2^M-j 



2nuj 



n 



Tb 



Tb 



Ui - U2 
Ui + U2 



- 1 



(144) 



where we have set ^—1, Ui = e^', Tb = e^^ cx A^/'^^^-'. 

Clearly, the integrals are now convergent at low energies, and we can send A to zero. Since the IR 
value of the one point function is easily determined by other means, < (7{x) >/_roc x^^^^ 



we can now 



obtain < a{x) >Tb from (144). Hence the procedure involves a double regularization. Of course, there 
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Figure 5: Accuracy of the finite Tb over the IR part of the envelope of p{x) for g = 1/2. 



remains an infinity of terms to sum over. However, as in the case of current operators, the convergence of 
the form-factors expansion is very quick, and the first few terms are sufficient to get excellent accuracy 
all the way from UV to IR. To illustrate this more precisely, we recall that for g = 1/2 (144) can be 
resummed in closed form p33, p2] , giving rise to : 



R 



'exact 



{(T[X))IR VTT 



(145) 



By reexponentianing the two first terms in ( |144|) , one gets a ratio differing from ( |145| ) by at most 1/100 
for xTb G [0,oo) (see figure 1). 

By reexponentiating the first three terms, acuracy is improved to more than 1/1000. Clearly, the form- 
factors approach thus provides analytical expressions that can be considered as exact for most reasonable 
purposes. 

It is fair to mention however that, at any given order in (144), the exponent controlling the a; — > 
behaviour is not exactly reproduced, as could be seen on a log-log plot. For instance, the first term is 
immediately found to produce a behaviour R{x) oc x^^^, to be compared with the result Rexact{x) oc 
x^/^lnx. The comparison of the exact result (145) and of (144) show that the form factors expression 



has, term by term, the correct asymptotic expansion ie the IR expansion in powers of jy^- Adding 
terms with more form-factors simply gives a more accurate determination of the coefficients. This is to 
be compared with the results of |35| for eg the frequency dependent conductance, where the form factors 
expression had the correct functional dependence both in the UV and in the IR. This is not to say that 
the method is inefficient in the UV, because we know, at least formally, all the terms. In fact, we will 
show in what follows how the expansion ( |l44|) can always be resummed in the UV, and that the exponent 
can be exactly obtained from the form-factors approach too. 

The regularization is the same for other values of g. Let us discuss here the cases g = 1/t with t 
integer again. For these values, the scattering is diagonal and the form factors are rather simple. To 
obtain them, we again take the massless limit of the results in |^ (but this time for vertex operator 
cos$/2 (in normalizations where the bulk sine-Gordon perturbation is cos$) instead of the current) 
and impose that half of the quasiparticles become right movers and half become left movers, since the 
boundary state always involve pairs of right and left moving particles. It is in fact easier to take that limit 
if we change basis from the solitons and anti-solitons to :^(|*S' > ±|A >). In that case, the boundary 
scattering matrix becomes diagonal and the isotopic indices always come in pairs. The reflection matrices 
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in this new basis are given by : 

R_{e) = -e'f (2-t) tanh + z^^j i?(if - 9) 

i?_^(6i) = e*f (2-t)i?(i| - 61) (146) 
with (note the shght change of notation compared with (|85|)) : 



R{e) = exp (^i J 



dy 2{t-l)y0 smhjt - 2)y 

2y TT sinh(2y)cosh(t- l)y/ ■ ^ ^ 



The breathers reflection matrices foflow from |Q (see also (p7|). 

The case g — 1/2 having aheady been worked out, let us concentrate on 17 = 1/3. There, in addition 
to the soliton and anti-soliton, there is also one breather. The first contribution to the one point function 
comes from the two breathers form-factor, with one right moving and one left moving breather. It is 
given by a constant : 

f{9,9)i,^^c,, (148) 

and this obviously leads to IR divergences. Other contributions come from 2n breathers form-factors, 
and 4n solitons form- factors The whole expression can be controlled as for g = 1/2, by taking the log, 
and factoring out the IR part. Setting c{x) — cos we organize the sum as follows : 

lni£Mlk^lni?(2)+lni?W + ... (149) 
{c[x))iR 

with the subscript denoting the number of intermediate excitations. 
Then, using the explicit expressions for g = 1/3 wc find : 

lni?(2) = 2cie^'^^'''' Ei{~2V2TBx), (150) 

where Ei is the standard exponential integral. The next term Ini?*^^' is a bit bulky to be written here, 
but it is very easy to obtain - similar expressions have been explicitely given in the previous sections. This 
is all what is needed for an accuracy better than 1 percent. In figure 6 we present the results of the ratio 
at gr = 1/2,1/3,1/4 for the Friedel oscillations. It should be noted that this ratio is just the pinning 
function of reference |2|] and our results agree well qualitatively with the results found there. 

As mentioned before, the deep UV behaviour is a little more difficult to obtain: the accuracy is good 
because the ratio goes to zero anyway, but the numerical evaluation of the power law would not be too 
accurate with the number of terms we consider. Fortunately, the full form-factors expansion allows the 



analytic determination of this exponent. First, observe for instance that in (144) the integrals converge 



for all Tb 7^ 0, but strictly at Tb = 0, they do not. To get the dependence of {c{x)) as Tb 0, we 
will consider, the logarithm of another ratio. In f^/^m^"" , where x and x' are two arbitrary coordinates. 
For this ratio, even at Tb = 0, the integrals are convergent. But Tb = is the UV fixed point, with 
Neumann boundary conditions. While the one point function (c(x))t/y vanishes, the ratio of two such 
one point functions is well defined, and can be computed by putting an IR cut-off (a finite system). One 
finds that it goes as {x/x'^^'^. By regularity as Tb 0, the same is true for the ratio close to Tb — 0, 
and thus one has 

(c(x)) oc {xTBy/^x{TB) ^ (151) 



This shows that the universal scaling function in (142) behaves as F{y) cx y^-^ for g < j. This exponent 



can actually be obtained by perturbation theory. Indeed, the first term in the perturbative expansion of 
{c{x)) is 

dy 

For g < 1/2, this integral diverges in the IR. To regulate it, we need to put a new cut-off: since there 
is no other length scale in the problem, this can be nothing but 1/Tg. Changing variables, the leading 
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behaviour is x^/'^Tg cx x^/^A^/^^^, in agreement with the previous discussion. The exponent coincides 
with the result of the self-consistent harmonic approximation ||2^ ; but it is important to stress that the 
latter is valid only for g << 1. 
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Figure 6: Ratio of the finite Tb over the IR part of the envelope of p{x).. 

The function F{y) behaves a j/lny for g — 1/2. For g > 1/2, its behaviour is simply F(y) cx ?/, as 
can be easily shown since the perturbative approach is now convergent. As we approach g — 1/2 this 
exponent seems to become asymptotic and is more difficult to get numerically p9| . 

7 Conclusion 

The description of the Hilbert space in terms of the integrable quasiparticles basis allows, in a wide 
variety of cases, extremely accurate computations of time and space dependent correlators. The related 
form-factors method is not especially elegant, but gives rise to susprisingly good results: in fact, it is not 
clear whether exact analytical expressions of the correlation functions - if ever obtained - will be more 
useful. Although we limited ourselves to g — 1/t with t integer in these notes, all values of g < 1/2 are 
accessible, but general computations are more complicated since the bulk scattering is non diagonal. The 
methods explained here should be generalizable to other problems, in particular the determination of the 
screening cloud in the anisotropic Kondo model . 

The region > 1/2 - in particular the SU{2) symmetric point g = 1 1^ presents additional difficulties, 
which are not yet solved. The situation does not look desperate however. 

As a last but important comment, I would like to stress that the form- factors technique has also been 
applied succesfuUy to the computation of quantities of experimental interest in the case of (bulk) massive 
theories. See the lectures by F. Esslcr in this school, or the papers [^, One point fimctions in 

massive theories and in the presence of a boundary have also been recently studied in B3| . 



Acknowledgments: These notes result from close collaborations with F. Lesage and S. Skorik, to whom 
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acknowledged. 

^Of course the problem of Friedel oscillations for g = 1 can be solved by fermionization [ pof , since it corresponds to non 
interacting fermions. In our approach, this point is non trivial because of the folding. This folding however is necessary for 
any value g 7^ 1: except at g = 1, the problem on the whole line would not be integrable. 
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